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ABSTRACT 



Laminar/turbulent transition is studied in a straight channel with 1 Hz 
imposed bulk flow unsteadiness for Reynolds numbers of 1 103 to 27 15 and 
Strouhal numbers of 0.06045 to 0.02455. Channel aspect ratio is 40 to 1 with 
1.27 cm height, 50.8 cm width and 4.27 m length. Observations, videos, and 
photographs of smoke patterns show different subcritical transition events 
including Tollmien-Schlicting waves, vortex-array type motion evidenced by 
smoke swirls and ribbon like patterns, turbulent spots and fully turbulent flow. 
Fluctuating intensity magnitudes relative to velocities from imposed 
unsteadiness are determined from phase-averaged velocity traces. Of 
particular interest, is whether the imposed unsteadiness advances or delays 
transition events. With the conditions that were investigated, transition 
events begin to occur at lower Reynolds numbers and extend over a wider 
range of Reynolds numbers compared to flows with no imposed unsteadiness. 
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L INTRODUCTION 



A. BACKGROUND 

Many typical flows found in nature and in aerodynamic applications and 
almost all flows associated with turbomachinery involve freestream 
unsteadiness. Most of these flows, especially those of technological interest 
undergo transition from a laminar state to a turbulent condition. The accurate 
prediction of the location and extent of transition and the accompanying 
changes of important flow properties is vital to the design of aircraft and most 
energy transformation machines involving fluid flow. 

Laminar-turbulent transition results in significant increases in local wall 
shear stress and convective heat transfer rates. This becomes evident by first 
considering laminar flows, where transport of heat and momentum normal to 
streamlines occurs only as a result of molecular diffusioa This is usually an 
advantage since tangential forces at solid surfaces (drag) and heat flux normal 
to surfaces are relatively low. Occasionally, laminar flows are undesirable since 
they tend to separate from surfaces more easily than turbulent flows. In most 
applications, laminar flows are unstable and transitional or fully turbulent flows 
exist In the latter, higher rates of diffusion across mean streamlines due to 
turbulent transport cause significantly higher drag and heat flux at surfaces. 

The parameter generally used to describe the onset of transition is the 
Reynolds number. In channel flows, linear stability theory predicts that the 
laminar flow becomes unstable to small amplitude disturbances at a critical 
Reynolds number of 5772. Here the Reynolds number is based on mean 



i 



centerline velocity and channel half-width. With large amplitude disturbances, 
however, the transitional Reynolds number for channel flow may be as low as 
1000- 1500. The value of Reynolds number at the onset of transition depends 
strongly on the disturbances existing in the flow and conditions which prevail 
in the entry region of the channel. 

Also, the sequence of events which occur during transition are strongly 
dependent on upstream history, especially laminar boundary layer character 
and development There is a growing collection of evidence which also 
suggests that turbulent shear layer development is more strongly dependent 
on its upstream history than previously believed. This would mean; ( 1) that 
transition affects turbulent boundary layer structure downstream, and (2) that 
turbulent boundary layer results obtained under equilibrium conditions (Le. 
streamwise self-similarity) may not adequately represent flow behavior near 
surfaces of operating turbines and airfoils. Non-equilibrium shear layer 
development near turbine passage surfaces and near airfoils is due to varying 
surface curvature, varying pressure gradients, and freestream unsteadiness. 
Unsteadiness is particularly important in this regard since it affects shear 
properties at a given location in two ways: (1) by its local direct influence and 
(2) from its affect on shear layers upstream, whose history affects downstream 
development The computational work of Singer, Ferziger and Reed [Ref. 1[ 
indicates that moderately low frequency oscillations in the channel can initiate 
nonlinear effects which trigger transition on flow regimes where the steady 
flow remains laminar. These happen even though linear theory indicates that 
the flow is stable. Experimental clarification of these points is needed, along 
with the verification of the sequence of events occurring during the transition 



process at different Strouhal numbers. The view of the overall transition 
process is especially valuable to the gas turbine research community. 

Unsteadiness is additionally important since boundary layers which exist 
near turbulent blades are significantly different from equilibrium boundary 
studied in many laboratories. A number of recent studies indicate the 
inadequacy of results obtained from equilibrium boundary layers. Gillis and 
Johnston's [Ref 2| measurements in boundary layers downstream of a convex 
curved surface and Ligrani's [Ref. 3) results in boundary layers downstream of 
artificial thickening devices showed that many structural changes persist long 
distances downstream of perturbations. This indicates that boundary layer 
behavior may depend on upstream history to a greater extent than previously 
believed This is important regarding unsteadiness since it alters boundary 
layers locally as well as the history of boundary layers upstream. 

Few experiments have been performed to study the stability or transition 
of oscillating flow in a plane channel. Most known information is from a 
handful of analytic studies and the full Navier-Stokes simulation undertaken at 
Stanford 

Regarding the latter, the work by Singer, Reed and Ferziger [Ref. 1] 
provides fascinating information as to the effects of unsteadiness on transition. 
For plane channel flow linear theory has been verified at a Reynolds number of 
5000. With mean flow frequency in the range of the non-dimensional 
frequency of the least stable TS wave (Tollmien-Schlichting), they find that 
the mechanisms driving the transition are similar to those described by 
Herbert's secondary instability theory [Ref. 4 and 5] and observed 
experimentally by Nishioka eL aL [Ref. 6 and 7], and Kozlow and Ramazanov 



[Ref. 8). In the plane channel flow these mechanisms have been studied 
extensively in previous computational investigations. The addition of mean 
flow modulation causes significant variations in the amplitude of the TS wave 
during the period, hence the dominant flow structures depend on the 
Strouhal number and the magnitude of the modulation. 

In a straight channel without imposed unsteadiness the process of 
transition begins with TS waves experiencing mild growth at sufficiently high 
Reynolds numbers. When these waves grow to an amplitude of about 0.3 
percent of the free stream velocity, three-dimensional waves become 
unstable. Transition often occurs at Reynolds numbers below that at which TS 
waves become unstable. These waves cause Lambda shaped vortices to lift 
from the near wall region and "burst" ejecting low speed fluid from close to 
the wall into the faster moving fluid away from the wall. These events then 
lead to a rapid breakdown of the laminar boundary layer into a turbulent one. 
The process can be modified by a number of different factors. Freestream 
turbulence can lead to non-sinusoidal oscillations in the nominally laminar 
boundary layer. Local disturbances, such as roughness elements, may lead to 
structures known as turbulent spots. These patches of turbulent-like flow 
form a characteristic arrowhead shape and spread in the downstream 
direction. Distributed roughness, acoustic disturbances, surface waviness and 
other factors all serve to complicate the basic process. There also exist 
equilibrium states where two dimensional waves (which are similar in form to 
TS waves) with a given amplitude will neither grow nor decay. The initial 
disturbances in the inlet of the channel can evolve into these or similar 
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structures, after which the strong secondary instability of three-dimensional 
waves can take place. 

B. OBJECTIVES 

This research is a follow-on to the work undertaken by Longest [Ref. 9], 
and is believed to be the first experimental (quantitative) study, known to the 
author, on the effects of imposed mean flow unsteadiness on transition in 
straight channel flow. The objective is to investigate the effects of imposed, 
periodic, and deterministic unsteadiness on the transition process from 
laminar to turbulent flow. Answers to the following questions are sought 

1.) What parameters govern the unsteady flow? 

2) What events occur in the flow during the onset and development of 
transition? 

3. ) What is the interaction between the organized induced unsteadiness 
and transition (e.g., development of two-dimensional Tollmien-Schlichting 
waves into three-dimensional instabilities)? 

4. ) At what range of frequencies can one expect interactions between 
the organized imposed unsteadiness and transition phenomena? 

Flow measurements (using hot-wire anemometry) and flow visualization 
are used to analyze the effect of oscillating the bulk flow in otherwise steady 
plane Poiseuille flow. An attempt is made to determine whether the type of 
unsteadiness studied (sinusoidal at one and two hertz) contributes to the early 
onset of transition, delay of transition, or has no effect. Sinusoidal oscillation 
was chosen to induce a deterministic periodic variation of velocity. The results 
are compared with the results of Singer's [Ref. 1] numerical simulation, which 



which suggests that such oscillations provide a stabilizing effect at all but very 
low frequencies. 

C. ORGANIZATION 

Section II follows this introduction and describes the experimental 
facilities. In section III experimental procedures are discussed, including 
calibrations, channel validatioa flow measurements and flow visualization. 

The phase averaging technique is also detailed here. A discussion of the 
results follows in section IV. Conclusions are presented in section V. 
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IL EXPERIMENTAL FACILITIES 



A. CHANNEL 

The 40 to 1 aspect ratio straight channel shown schematically in Figure 
1, is located in the laboratories of the Department of Mechanical Engineering. 
The details of the design and construction are discussed by Ligrani and 
Subramanian [Ref 10). Briefly, the channel is made of 6.35 mm thick 
plexiglass with a straight test section 4.27 m in length, with inside 
dimensions of 1.27 cm in height and 50.8 cm in width. The top and bottom 
walls of the channel are supported by ribs and cross beams along its length. 
The side walls are removable in order to gain access to the inside of the 
channeL The inlet section consists of a honeycomb assembly, 3 frames with 
screens and a 20: 1 contraction ratio nozzle. Two layers of cheesecloth are 
placed at the inlet lip for filtering the air entering the channel At the exit of 
the 4.27 m long test section are three 10.16 cm long frames with screens for 
housing the unsteady device ( the device used to create unsteadiness in the 
flow) and a second honeycomb. A two dimensional diffuser 45.72 cm long 
with a 3 degree total angle is located just upstream of the plenum chamber. 
The exit plenum chamber has inside dimensions of 60.96 cm by 60.96 cm A 
3.81 cm orifice plate assembly between the plenum chamber and the blower is 
provided to meter the flow through the channeL The exit of the plenum 
chamber is connected to the suction side of a 5 H.P. blower with a long 50.8 
cm diameter pipe. A valve is provided between the plenum and the blower to 
regulate the flow. Bulk flow velocities up to 15 m/s are achieved in the 
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channel. In our experiments, bulk flow velocities in the range 1.0-4.5 m/s 
are used 

The channel is designed with a sufficient length such that transition 
occurs after the laminar velocity profile (parabolic) is fully developed (self 
similar). The initial flow development length is estimated [Ref 1 l:p. 186] to 
be within about 25 per cent of the test section length for the range of speeds 
used in the present tests. 

B. UNSTEADY DEVICE 

The most important aspect in the design of any unsteady device for a 
transition study is to ensure that the device produces a deterministic and 
periodic unsteadiness without adding other disturbances. To achieve this in 
an open circuit induction tunnel the unsteady device is best located at the 
exit of the channeL This way, problems caused from the wakes and the flow 
disturbances resulting when the unsteady device is upstream of the test 
section are avoided. 

The design is based on the design used by Miller and Fejer [Ref. 12] and is 
shown in Figure 2 Since the depth of the channel is only 127 mm, a single 
rotating vane driven through a spur gear train by a DC stepper motor is 
adequate to introduce the required unsteadiness in the flow. The rotating 
vane causes varying resistances to the flow in the test section. The frequency 
of the imposed unsteadiness is controlled by changing the vane rpm and the 
amplitude is controlled by using different vane widths. It is also important to 
match the flow speed with the vane speed. A vane 8.7 mm wide was used in 
this study. 
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The vane in Figure 2 is made of a 3.2 mm thick brass strip with rounded 
edges, and spans the entire width of the channel. The vane is supported at 
the ends by a 3.2 mm diameter shaft and bushings that are fitted to the side 
walls of the frame. Three additional intermediate spanwise struts are 
provided to increase the rigidity of the vane. One end of the shaft is extended 
to accommodate a 48 TPI (Threads Per Inch) spur gear. The driver spur gear 
has 12 TPI and is mounted on the Superior Electric, M092-FD310 Stepper 
Motor. To achieve a fine step control, a Modulynx MITAS PMS085-D050 
Drive with Modulynx MITAS PMS085-C2AR Drive Controller is used. This 
combination facilitates motor shaft speeds up to 100 revolutions per second, 
with a constant 200 steps per motor shaft revolution. Note that, with a 1:4 
gearing ratio and two cycles of imposed flow oscillations per vane rotation, one 
motor revolution corresponds to one half cycle of the imposed flow oscillation 
by the vane. The Drive Controller can be programmed for the required motor 
speed (rps). The Modulynx PMS085-D050 Drive was also used to provide a 
trigger signal to initiate data acquisition. This trigger also enabled the 
correlation of data taken with the position of the vane at the time that the data 
was recorded 

C. FLOW MEASUREMENTS 

The equipment used for flow measurement includes a digital manometer, 
hot-wire probe, hot-wire bridge, signal conditioner, digital multimeter, high 
speed data acquisition system, oscilloscope and a microcomputer. These are 
described in detail below. 
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L Bulk Velocity 

The bulk flow velocity through the channel is determined from the 
pressure drop across the 3.8 1 cm orifice plate which is installed in the piping 
between the blower and the exit plenum chamber. This pressure drop is 
measured and displayed by a Validyne Model PS 309 digital manometer 
(inches of water). 

Z Instantaneous Velocity 

A constant temperature hot-wire anemometer is used to measure 
instantaneous velocities of the flow. Briefly, the principles of hot-wire 
anemometrv (constant temperature mode) are as follows. 

The flow sensing element (hot-wire probe) consists of a short length of 
wire, connected to two supports, which forms the leg of a Wheatstone bridge. 
This bridge maintains the wire at elevated temperature and is initially 
balanced prior to exposing the probe to the flow. This balancing resistance, 
referred to as "cold" resistance is multiplied by an "overheat ratio". Once the 
bridge is set to this value, and as the probe is exposed to the flow, the bridge 
is placed in operation and a small current (60 mA) passes through the circuit 
Wire temperature increases because of resistive heating and the wire 
resistance increases until the bridge is rebalanced. The governing relationship 
is: 

R = Ro(l +a(7 - To)) 

where R is the resistance at temperature T, Ro is the resistance at the 
reference temperature To, and a is the temperature coefficient of resistance. 



Flow past the probe causes the wire to cool (by convection) and the 
bridge becomes unbalanced as the resistance decreases. The bridge is 
rebalanced by a feedback loop and amplifier. Response is rapid enough to 
maintain wire temperature and resistance nearly constant As such, the 
variations in the bridge voltage are directly proportional to the current 
variations, and flow velocity can be calculated as a function of the bridge 
voltage. 

Now the different components used to measure instantaneous velocity 
are described. 

a. Hot-wire Probe and Probe Positioning 

A Dantec Electronics Inc. P51 probe, mounted such that the wire is 
horizontal and normal to the streamwise flow direction, is used. The wire 
sensor is 5.0 microns in diameter and approximately 1.25 mm long. The 
probe is mounted through the side wall of the channel and is held fixed by a 
rotatable lever arm which allows placing of the probe at various y/d positions 
ranging from 0.2 to 0.9. A scale graduated in terms of y/d is mounted adjacent 
to the lever arm for probe positioning to an accuracy of approximately 0.5 mm 
Longitudinal and transverse position of the probe for this study were 46.12 cm 
and 15.64 cm, measured from the near flange at the end of the test section, 
and from the inside of the left channel wall (looking downstream), 
respectively. By rearranging the sidewalls the probe can be positioned at 
other streamwise locations for future studies. 

b Hot-wire Bridge 

The hot-wire is operated with a DISA 55 M 10 constant temperature 
bridge. A five meter coaxial cable is used to connect the hot-wire to the 






bridge. For the purposes of this study the cold resistance was multiplied by an 
overheat ratio of 1.8. The DC voltage from the bridge is measured with a 
Hewlett-Packard 3466A digital multimeter. 
c. Signal Conditioner 

The hot-wire signal is filtered and amplified by a DISA 56W20 signal 
conditioner. The high-pass filter is set to 0.1 Hz in all cases to remove the DC 
signal This however, did not alter the effects of 1 and 2 Hz imposed flow 
unsteadiness. The low-pass filter is set at 10 KHz to filter out noise existing 
above this frequency. The hot-wire signal is amplified using a gain of 10 to 
achieve an acceptable signal to noise ratio. Output from the signal conditioner 
is sent to the high speed data acquisition system and to a B & K Precision 10 
MHz Oscilloscope used to signal monitor the signal 
d High Speed Data Acquisition System 

A Hewlett-Packard 6944A Multiprogrammer is used in conjunction 
with a Hewlett-Packard series 9000, Model 310 microcomputer to acquire 
data. A buffered 69759A A/D card configuration is used with the 
Multiprogrammer to allow continuous data acquisition. The A/D card is 
externally triggered using 2-3 volt TTL type signals from the unsteady device 
motor drive. When the unsteady device is not in operatioa triggering is 
provided by a Hewlett-Packard 33 1 LA Function Generator set to the same 
frequency as that used with the unsteady device in operatioa Sampling 
occurs as square wave voltage changes from high to low. 

The storage capacity of the A/D buffer is 64,000 bytes. The number of 
imposed unsteadiness cycles which can be sampled depends on the desired 
number of data points per cycle. For example, if 100 points per cycle are 



desired, then 640 cycles can be sampled. In this study, sampling was done at 
100 points per cycle for all cases, where cycle refers to the cycle of flow 
blockage caused by the vane rotation of the unsteady device. 100 points per 
cycle were chosen so that there would be sufficient resolution to determine 
the effect of unsteadiness over the complete phase of each unsteadiness cycle. 
e. Datastorage 

The Hewlett-Packard microcomputer is used to acquire the data from 
the A/D buffer, using a data acquisition program HOTWIREPAV, described in 
the software directory (APPENDIX A). Raw instantaneous voltage data is 
acquired and stored on standard 3.5 inch micro diskettes for further 
processing. Manual inputs to the program include bulk flow velocity, flow 
blockage frequency, and date and time. 

£ Data Processing 

Processing of the data is done using a Hewlett-Packard 
microcomputer using the programs MAX_M1N and AVGPLOT. This software is 
described in detail in the software directory (APPENDIX A). Many of the 
curves used to graphically display the results were generated using a 
Macintosh SE personal computer and the graphics program CRICKETGRAPH. 
The phase averaging technique employed by the program MAX_MIN is 
described below. 

In the case of an unsteady flow the instantaneous velocity can be split 
into three components. 

U » D + U + o 

where V is the time averaged velocity, u is the imposed phase averaged 
velocity and if is the instantaneous fluctuating component u is the 



component due to the unsteadiness, and is zero for steady bulk flow. In the 
present study u and V are combined to represent the phase averaged velocity 
{ (J). Thus, 

U = u + (/ . 

u is determined by phase averaging U over, say 640 cycles (ncycl), at 
different positions of the phase of each oscillation. 
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When if is a large percentage of u then a larger number of cycles, 
sometimes up to 4640 cycles, are used to determine U . For this, multi-run 
data acquisition is made, with the vane at the same starting position for each 
run. Only the last 580 data points (due to sampling procedure) of each run are 
used in the cumulative phase averaging. Once h is determined, if is 
calculated from if - U - u , and ((/)" - (U - i /f . Therefore, the phase 
averaged root mean squared velocity is, 
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Time averaged velocity, 0 1 and rms velocity are calculated by 

averaging over all the points (100 x ncycl) per run and in the case of multi- 
runs, 100 x ncycl x no. of multi-runs. 



\A 



«wnV 



z?« 



y 

ncvcf 

/n-1 



- IW 

ioo^.") 



/> -1 




J V 



00 



/?) - r/)‘ 

/>- 1 



D. FLOW VISUALIZATION 

The equipment utilized for flow visualization studies includes a smoke 
wire for generating smoke in the channel and video equipment to photograph 
the flow in the channel. These are described in detail below. 

1. Smoke Wire System 

The smoke wire system is a modified Vertical Smoke Wire Instrument 
made by Flow Visualization Systems of Bolingbrook, IL. The system in its 
present configuration consists of a 0.127 mm diameter nichrome wire that is 
stretched across the width of the channel through access ports in the channel 
side walls. The wire is passed through these access ports through a thin 
stainless steel tube (hypodermic needle) which allows the wire to move freely. 
The smoke wire is connected to an A/C-D/C converter and this in turn has 
voltage supplied by a Calrad 45-740 (0-130 VAC) variac set at 45 VAC. The 
wire is energized by pulsing the ’Momentary On" switch of the A/C-D/C 
converter. Since the wire acts as a resistor and heats up when a current is 
passed it becomes more elastic and has a tendency to sag in the channel To 
prevent this from affecting the flow visualization, a tensioning device has been 
installed on one end of the wire to ensure that the wire is held taut in the 
channel when it is energized. This tension device consists of two pulleys 






which support the wire and a third pulley with a weight attached between the 
two support pulleys. The wire is coated with a paraffin based oil which is 
applied by means of a thin rod with an absorbent piece of felt at the tip. When 
applied in this manner, the wire is coated evenly with small droplets of oiL 
Energizing the wire produces a thin even sheet of smoke until all the oil on 
the wire bums. The smoke from the wire lasts for about 4-8 seconds 
depending on the flow rate in the channel. 

The wire is placed horizontally normal to the flow 76.8 cm upstream 
from the end of the test section. For the purposes of this study, the wire was 
located at a y/d value of 0.84 (d being the width of the channel and y the 
height from the bottom of the channel). 

2. Video Equipment 

The video equipment used for flow visualization is the same as that 
described by Longest (Ref. 9). Briefly, a high speed video camera is mounted 
on a tripod above the channel and focused on the area of interest The focal 
plane from the camera lens to the channel is surrounded by an enclosure to 
exclude all extraneous light and thus decreasing the reflections which 
interfere with the flow visualization. The bottom of the channel was blacked 
out and lighting was provided by two 600 Watt floodlights placed 1 m from the 
sidewall on either side of the channel at the same height as the channel 
Still photographs were made by taking pictures of the aforementioned 
video when viewed on a monitor. A Nikon F-3 SLR and KodalcTri-Xpan film 
were used to take the pictures. The shutter speed was one-eighth of a second 
and the focal aperture was set at f 5.6. 



ILL EXPERIMENTAL PROCEDURES 



A. ORIFICE AP CALIBRATION 

The orifice differential pressure calibration for bulk velocity flow 
measurements was performed by Longest [Ref. 9). Results of the calibration are 
shown in Figure 4. 

R HOT-WIRE CALIBRATION 

The hot-wire was calibrated in the freestream of the wind tunnel, also 
located in the laboratories of the Mechanical Engineering Department of the 
Naval Postgraduate School. 

The hot-wire probe is mounted normal to the flow in the center of the 
channel so that it is in the freestream. The hot-wire is connected to the same 
hot-wire bridge previously discussed. Output from the hot-wire bridge is read 
on a digital multimeter. Freestream velocity in the wind tunnel is measured 
utilizing a Kiel pressure probe, a wall static pressure tap, and a Validyne digital 
manometer. Voltage and differential pressure readings are taken for a range 
of pressure drops corresponding to velocities of 1.0 to 4.0 m/s. The computer 
program HWCAL is run on the Naval Postgraduate School's IBM main frame 
computer to calculate all calibration constants. Following the calibratioa the 
hot-wire is installed in the straight channel. 



C. CHANNEL VALIDATION 



1. Channel Flow Tests 

Assembly and installation of the straight channel in the Mechanical 
Engineering Laboratory of the Naval Postgraduate School, along with the initial 
flow visualization and measurement results are detailed by Longest |Ref. 9]. 
Prior to commencing flow measurements for this study a check of the channel 
flow quality is performed as follows. Smoke is generated using the smoke 
generator described by Longest [Ref. 9] and then drawn into the channel at a 
speed of approximately 20 m/s. The first few times this was performed, 
numerous leaks from the sidewalls of the channel were observed These leaks 
are eliminated using removable spring clamps on both sides of the channel, 
along its entire length. The clamps are left in position while the channel is in 
operation. 

2 Laminar Average Velocity Profile Verification 

Velocity profile verification is performed at a Reynolds number of 1103 
(bulk velocity of 1.32 m/s) when the flow is laminar. Instantaneous velocity 
measurements were taken at eight y/d positions ranging from 0.2 to 0.9, at an 
x/d of 300.0. This data was taken using the HOTWIREPAV program, with no 
imposed unsteadiness. MAX_MIN was used to process the data, and the time 
averaged local velocities were normalized with respect to the maximum value 
which occurred at the mid-channel location (y/d = 0.5). A second order 
polynomial curve was fit to the data and is shown in Figure 5. The curve 
demonstrates that the flow is fully developed at this location and follows the 
expected parabolic behavior. 
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D. FLOW MEASUREMENT 
1. Reynolds Number Survey 

Velocity surveys are conducted for Reynolds numbers ranging from 1 100 
to 2700 with the hot-wire probe at a y/d position of 0.85, and a z/d of -7.7. 
This survey is done to determine how the imposed unsteadiness affects the 
transition process as Reynolds number varies. Data at each Reynolds number 
are taken by recording results with no imposed unsteadiness followed by 
results with imposed unsteadiness. Data sets of 64,000 points are taken using 
only the last 58,000 for data processing, as previously described. Average 
velocity and rms velocity values from such survey are used to generate Figures 
6 and 7. These experiments were repeated at several Reynolds numbers to 
ensure consistent and repeatable data trends. 

Z Multi-run Averaging 

Multi-run phase averaging using up to 4640 cycles of flow unsteadiness 
were conducted at several Reynolds numbers. The purpose of this multi-run 
averaging is to show that the survey data is statistically stable for higher 
Reynolds numbers where fluctuation levels are significant with respect to the 
local mean velocity. Results are presented in Figures 8-35. 

3. Velocity Profiles 

Velocity profiles with and without imposed unsteadiness are taken at the 
same Reynolds numbers as used for the multi-run averaging. The purpose of 
these profiles is to record flow development as the Reynolds number is 
changed and to study the effect of imposed unsteadiness on average and rms 
velocities as a function of the Reynolds number and y/d. Data acquisition and 



processing is performed as described in the Reynolds number survey. Results 
are provided in Figures 36-147. 

E. FLOW VISUALIZATION 

Smoke patterns were recorded on film using the video camera at Reynolds 
numbers ranging fromll40 to 2240 with 1 Hz imposed unsteadiness. Video 
sequences are taken at different Reynolds numbers spaced apart by 
approximately 50 units. The purpose of this flow visualization is to record 
transition events and denote any variations caused by the imposed 
unsteadiness. Still photographs are shown in Figures 148-158. 
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IV. RESULTS AND DISCUSSION 



The coordinate system used in this study is a right handed coordinate 
system. The origin is at the spanwise centerline, at the bottom of the channel 
and at the beginning of the test section. The x-axis extends positively in the 
streamwise direction to the end of the test section. The y-axis extends 
vertically in the positive direction from the bottom of the channel to the top. 
The z-axis extends in the spanwise direction. Positive values of z are on the 
right hand side of the channel when looking downstream. Flow 
measurements were taken at x/d = 300, zfd = -7.7 and y/d values of 0.2 to 0.9, 
where d is the full channel height Flow visualization was done with the 
smoke wire positioned at a y/d value of 0.85. 

The Reynolds number. Re. is based on the bulk flow velocity, Vbujk, and the 
channel height and is defined by the equation: 

Re = Vtour • d/v 

where v is the kinematic viscosity. The Strouhal number, Str, is defined by 
the equation: 

Str = 2a* fcec* d/vfcuik 

where fcec is the frequency of the imposed unsteadiness. This study was 
conducted at Reynolds numbers from 1103 to 2715. With 1 Hz imposed 
unsteadiness, Strouhal numbers range from 0.0605 to 0.0246. 



A. CHANNEL FLOW QUALIFICATION 

For fully developed channel flow the theoretical velocity relationship is 
given by Schlichting [Ref ll:p. 280-281] as: 
u/umwt = 4.0[ y/d - (y/d) 2 ] 

where u is the local velocity and Udm* is the centerline (y/d = 0.5) velocity. 
Figure 5 shows the velocity profile taken at x/d = 300 at a Re of 1 103. The 
velocity profile (normalized with respect to Umax) is parabolic in shape and is 
in good agreement with the theoretical relationship. With the bulk flow 
velocity of 1.32 m/s (measured from the orifice pressure drop) the maximum 
centerline velocity of 1.774 m/s is observed, giving a ratio of Vb^/umax of 
0.774. 

B. EFFECT OF REYNOLDS NUMBER ON LOCAL MEAN VELOCITY AT 
y/d = 0.85 DUE TO 1 HZ IMPOSED UNSTEADINESS 

Figure 6 shows the variation of the mean velocity with Reynolds number 
with and without imposed unsteadiness. Flow measurements are for y/d = 
0.85. Bulk velocity is also shown in the plot. The mean velocity trace shows 
some non-linear behavior for Reynolds numbers from 1900 to 2200 and near a 
Reynolds number of 2700. This behavior is believed to be caused by the 
transitional behavior of the flow. The 1 Hz imposed unsteadiness has no effect 
on the mean velocity. 



C. EFFECT OF REYNOLDS NUMBER ON RMS VELOCITY AT y/d = 0.85 

Figure 7 shows normalized longitudinal velocity fluctuations with and 
without imposed unsteadiness for Reynolds numbers from 1 100 to 2700. 

The time averaged rms velocity, -yfu^ , is normalized with respect to the time 
averaged velocity, It , and is plotted versus the Reynolds number. For 
Reynolds numbers of 1100 to 1900, the imposed unsteadiness appears to have 
no effect on the rms velocity. For Reynolds numbers from 1900 to 2025, the 
1 Hz imposed unsteadiness causes a significant increase in the rms velocity. 
This increase amounts to a maximum at a Re = 1980 and is 7.06 %. At 
Reynolds numbers of 2025 to 2300 the effect of the imposed unsteadiness is 
reversed. In this range the imposed unsteadiness reduces the normalized rms 
velocity. The maximum amount of decrease occurs at Re = 2172 and 
corresponds to a drop of 7.3 %. At higher Reynolds numbers, the imposed 
unsteadiness appears to have no effect on the rms velocity. The figure also 
seems to demonstrate that transition begins near Re = 1200, and that the flow 
becomes fully turbulent at Re = 2400. 

D. EFFECT OF REYNOLDS NUMBER ON PHASE AVERAGED VELOCITY 

Figure 8 through Figure 35 show the phase averaged velocity traces at five 
Reynolds numbers (1153,1826, 2005, 2085, and 2640). Figures 8-17 are the 
phase averaged plots with no imposed unsteadiness. Of these, figures ( 8. 10, 

A, A 

12, 14, and 16) plot phase averaged velocity in this case U = D since U = 

0) versus phase angle. Figures 9, 11, 13, 15, and 17 plot the phase averaged 
rms velocity, V u ' ~ , normalized with respect to the temporal mean, D , versus 
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phase angle. For these cases, phase averaging is done over 580 cycles of 
imposed unsteadiness with 100 data points per cycle. 

Figures 18-35 present results of flow measurements with 1 Hz imposed 
unsteadiness. These figures show that the unsteady device produces 
deterministic and periodic unsteadiness, thereby validating its behavior. At Re 
= 1 153 the phase averaging is given for 580 cycles of imposed unsteadiness. 
For Re = 1826, 2005 and 2085 eight runs of 580 cycles each (for a total of 
4640 cycles) are averaged. For Re = 2640 three runs of 580 cycles (total of 
1740 cycles) are averaged. With 1 Hz imposed unsteadiness, statistically 
stable phase averaging is achieved using 580 cycles when the flow is laminar 
(Re = 1153 for the plots presented). However, for the transitional Reynolds 
numbers. Re = 1826, 2005 and 2085, eight sets of 580 cycles are necessary to 
achieve statistically stable phase-averaged results. For the turbulent case, Re = 
2640, statistically stable results are achieved after phase averaging three sets 
of 580 cycles of imposed unsteadiness. The effect of multi-run averaging is 
shown in Figures 20-35. Figure 20 shows U after averaging using one set of 
580 cycles. Figure 22 shows U after averaging using eight sets of 580 cycles, 
and is clearly much smoother indicating its statistical stability. Figures 21 and 
23 show the same effect on the rms velocity. Figures 24-27 show the same 
results of multi-run averaging for Re = 2005, Figures 28-31 for Re = 2085 and 
Figures 32-35 for Re = 2640. 

Normalized rms velocity values in these figures are fairly constant as the 
phase angle changes at all Reynolds numbers, with the exception of Re = 

2085. These data are shown in Figure 31. Here, the rms velocity increases 
for the phase angles of 250-350. This increase occurs at the same phase 



angles where u is larger than surrounding values. The increase in rms 
velocity values may occur during decelerating portion of the h cycle, followed 
by a phase shift relative to the u trace. These rms velocity variations occur 
over the same Reynolds number range that imposed unsteadiness decreases 
the rms velocity relative to the flow without imposed unsteadiness. 

EL MEAN VELOCITY PROFILES AND PHASE AVERAGED VELOCITIES 

Figures 36-42 show mean velocity profiles at x/d = 300 and z/d = -7.7. 
y/d positions vary from 0.2 to 0.9. Data is not given for y/d of 0. 1 due to 
physical constraints imposed by the geometry of the hot-wire probe. 

Figure 36 shows the profile development with Reynolds number with no 
imposed unsteadiness. Profiles are shown for Re = 1153, 1826, 2026, 2085 
and 2640. With laminar flow, profiles are parabolic. Profiles are somewhat 
contorted at transitional Reynolds numbers (Re = 1826, 2026 and 2085). At 
Re = 2640, the profile shows a flattened behavior expected for fully turbulent 
flow. 

Figure 37 shows the profile development for the same Reynolds number 
with 1 Hz imposed unsteadiness. The imposed unsteadiness nominally has no 
effect on the mean velocity profiles. This is also shown in Figures 38-42 
where mean velocity normalized with respect to the maximum observed mean 
velocity with and without unsteadiness are plotted for each Reynolds number, 
both with and without imposed unsteadiness. 

Figure 43 shows the variation of the peak to peak magnitude of the phase 
averaged mean velocity as a function of y/d. On the plot, the peak to peak 
magnitude has been normalized with respect to V . From the figure, the peak 
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to peak magnitude increases when the flow is near the wall At the two 
transitional Reynolds numbers (Re = 2026 and 2085), increases in magnitude 
near the wall are significantly higher than for other cases. These Reynolds 
numbers correspond to the situation in which the imposed unsteadiness 
decreases the rms velocity relative to flow with no imposed unsteadiness. 

Figures 44-9 1 show phase averaged data corresponding to data just 
discussed. These data indicate that the unsteady device is successful in 
producing deterministic periodic unsteadiness at all positions where flow is 
measured, i.e. from y/d = 0.2 to 0.9. 



F. RMS VELOCITY PROFILES 

Figures 92-98 show the development of the rms velocity profiles as the 
Reynolds number. Again, the symmetry of the profiles is distorted due to the 
fact that no measurements are available at y/d = 0.1. 

Figure 92 presents rms velocity profiles for the five Reynolds numbers 
with no imposed unsteadiness. Figure 93 shows the profile development with 
imposed unsteadiness at 1 Hz. Figures 94-98 are normalized rms velocity 
profile plots with and without imposed unsteadiness for each Reynolds 
number. 

Figures 92 and 93 show that the imposed unsteadiness has no effect on 
the rms velocity profiles for the laminar and turbulent Reynolds numbers (Re 
= 1153, 1826 and 2640). For the turbulent case, the rms velocities are much 
higher near the wall For the two transitional Reynolds numbers, some 
different behavior is seen. With no imposed unsteadiness, the profiles show 
higher rms velocity magnitudes near the middle of the channeL With imposed 
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unsteadiness, the profiles at Re = 2026 and Re = 2085 differ sharply and 
resemble the turbulent profile. 

Figures 99-146 present phase averaged data used to produce the plots 
described above. Rms values are constant as a function of phase angle both 
with and without imposed unsteadiness except for the case with imposed 
unsteadiness at Re = 2085. In Figures 116, 118, 120, 122, 124, 126, 128, 
and 130 rms velocities vary with phase angle. This was mentioned earlier for 
y/d = 0.85, however, here it is evident that this behavior occurs for 0.2 < y/d < 
03 . Figure 147 shows the variation of rms peak to peak magnitude as a 
function of y/d for Re = 2085. At this Reynolds number, the peak to peak 
magnitude is largest near the wall with a secondary peak near the channel 
centerline. 

G. FLOW VISUALIZATION 

Figures 148-158 show flow visualization results obtained with the smoke 
wire at y/d of 0.85. The camera was set up to take videos on the left side of 
the channel in the same area that the hot-wire was located (x/d = 300). In 
the figures, flow is from the bottom of the picture to the top. The left side of 
the picture is near the left channel wall, and the right hand side is located 
near the center of the channel. The dark horizontal stripe emanating from 
the left sides of the pictures indicates the location of the hot-wire probe. 

There are also several dark shadows which extend diagonally upward from the 
left side of the channel. These are due to shadows cast by part of the support 
structure of the channel and should not be confused with effects of the flow on 



>7 



the smoke patterns . The smoke generated by the wire appears as white lines 
or ribbons and the dark areas indicate absence of smoke. 

Figure 148 shows flow at Re = 1 136 with 1 Hz imposed unsteadiness. In 
this picture, the flow is fully laminar as shown by the straight regular parallel 
smoke lines. Figure 149 is at Re = 1 153 also with imposed unsteadiness. 

Here flow is again laminar, but there are also TS (Tollmien-Schlichting) waves 
in addition to some small streamwise spirals due to secondary vortical flows. 
These are hereafter referred to as ribbons. Figure 150 at Re = 1278 again 
shows laminar flow with TS waves. Figure 151 at Re = 1328 also shows 
laminar flow with TS waves. Figures 152 and 153 are at Re = 1395, and show 
TS waves and ribbons which are larger than previously seen and which begin 
to spread out diagonally. Figures 154 and 155 are at Re = 1428 again show 
the same ribbon structure, but here they have become larger and seem to 
develop into turbulent spots. These spots are defined by the large circular 
areas in the center of the ribbon break-up regions. Figures 156 and 157 are 
for Re = 1586. These pictures show the decay of ribbon structures as well as 
larger and more frequent turbulent spots. Figure 158 shows a time sequence 
of six pictures taken approximately one second apart for Re = 1428. In this 
sequence we see laminar flow w’ith TS waves is present in addition to ribbons, 
ribbons developing into turbulent spots, followed by more ribbons and laminar 
flow with TS waves. The principal difference between the flows with and 
without imposed unsteadiness is the reappearance of such TS waves at higher 
Reynolds numbers with imposed unsteadiness. In this situation, the TS waves 
generally reappear just downstream of turbulent spots. 



V. CONCLUSIONS 



The 1 Hz imposed bulk flow unsteadiness has no effect on the local mean 
velocity at all y/d values for the range of Reynolds numbers 1 100 to 2700. 

The imposed unsteadiness causes both destabilization and stabilization of 
the time averaged longitudinal velocity fluctuations. The flow is destabilized by 
the imposed unsteadiness for Reynolds numbers 1900 to 2025 and Strouhal 
numbers from 0.0351 to 0.0329. Stabilization of the flow by the unsteady 
device occurs for the Reynolds numbers 2025 to 2300 and for Strouhal 
numbers from 0.0329 to 0.0290. 

Imposed unsteadiness affects the shape of the time averaged longitudinal 
velocity profile at the transitional Reynolds numbers, Re = 2026 and Re = 

2085. 

The imposed unsteadiness also causes the time averaged longitudinal 
velocity fluctuations to vary as a function of phase angle at a Reynolds number 
of 2085 at y/d positions of 0.2 to 0.9. The maxima of the time averaged 
longitudinal velocity fluctuations and the maxima of the phase averaged 
velocity occur for the same phase angles. 

Transitional flow structures which have been predicted by numerical 
simulations and observed in other experimental work are also observed in this 
study. As Reynolds number increases laminar flow, Tollmien-Schlichting 
waves, vortical ribbon-like structures and then turbulent spots are observed. 
These structures are observed to be present at many different Reynolds 



numbers. Time sequence photography taken at particular Reynolds numbers 
during transition shows TS waves, vortical ribbons, turbulent spots, more 
vortical ribbons and TS waves. The growth and break-up of the ribbon-like 
vortical structures is followed by turbulent spots. As the Reynolds number is 
increased, the vortical structures grow in both the spanwise and streamwise 
direction. The frequency and size of the turbulent spots also increases as 
Reynolds number increases until the flow becomes fully turbulent The 
principal difference between the flows with and without imposed 
unsteadiness is the reappearance of the TS waves at higher Reynolds numbers 
with imposed unsteadiness. In this situation, the TS waves generally reappear 
just downstream of turbulent spots. 



APPENDIX A 



SOFTWARE DIRECTORY 

1. HWCAL: This program determines the constants for the King's Law 
calibration of the hot-wire. The program also provides a polynomial fit of the 
calibration data 

2 HOTWLKEPAV : This program is used to read the data stored in the A/D 
buffer of the high speed data acquisition system and stores the information on 
micro diskettes. Manual inputs are: triggering frequency, hot-wire DC voltage 
(ungained), oscillation frequency (flow blockage), bulk velocity and date and 
time of run. 

3. MAX-MIN : This program calculates instantaneous and phase averaged 
velocities. Initially a look up table is created. Here, effective velocities are 
calculated from the effective voltage values and stored for follow on 
calculations. The hot-wire calibration constants obtained from HWCAL, and 
the amplifier gain are incorporated into these calculations. The velocity 
calibration is given by the equation: 

UefT= k(Eeff2-Eo2)l/N 

where k is the proportionality constant, E e ffis the effective voltage, and E 0 is 
the reference voltage at no flow. N is a constant value of 0.45 for moderate 
Reynolds numbers. Once the look up table is created the program reads the 
instantaneous voltage values from the data file and converts them to 



instantaneous velocities. At this point a plot of the instantaneous velocity 
versus time can be generated for any of 580 cycles (The first 60 cycles are 
discarded to allow for flow stabilization as the unsteady device is started) 

Next the program phase averages the 580 cycles, and velocity versus phase 
angle (of the flow blockage) plots are available. Two plots are available from 
the averaged values; u versus phase angle and divided by U, versus 

phase angle, where u is the phase averaged velocity, VF is the phase 
averaged root mean squared velocity, and U is the average velocity. In the case 
where there is no imposed unsteadiness, the phase averaged velocity, u, is 
equal to the time averaged velocity , u, since there is no superimposed phase 
average velocity component, u. 

The program can also be used to process a set of data when multi-run 
averaging is desired In this case the averaged values are written to and stored 
in the computer hard disk memory. Once each of the desired data sets has 
been averaged and stored in the computer memory. 

4. AVGPLOT : The program AVGPLOT is used to average any number of data 
sets. A maximum of eight sets of data (or 4,640 cycles) were averaged for this 
study. The program AVGPLOT provides the same graphs as those available 
from MAX_MIN. Phase averaged data sets are copied to micro diskettes prior 
to purging them from the computers memory. 
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APPENDIX B 



FIGURES 



Figure 1 
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Figure 2. 
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Figure 3. 
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Figure 4. 
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Figure 11. 
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Figure 12. 
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Figure 31. 
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Figure 33. 
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Figure 34. 
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Figure 35. 
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Figure 37. 
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Mean Velocity Profile at Re No. = 1 153 
with 1 Hz Imposed Unsteadiness 
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Figure 38. 
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Mean Velocity Profile at Re No. = 1826 
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Figure 39 . 
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Mean Velocity Profile at Re No. = 2026 
with 1 Hz Imposed Unsteadiness 
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Figure 40. 
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Figure 41, 
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RMS Velocity Profile Development 
with No Imposed Unsteadiness 



o 


o 


O 


o 


o 


> 


> 


> 


> 




bl 


w 


w 


bl 


bl 


s 


s 


s 


6 


s 


to 




to 


to 


to 


00 


CD 


CD 


If; 


O 


If; 


CO 


CO 


CO 






CD 


o 


o 


CD 






CO 


CO 


CO 


■ 


□ 


< 


<1 


• 




*3 



Figure 92 . 



(s/ra) pisuu 



125 



RMS Velocity Profile Development 
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Figure 93. 
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RMS Velocity Profile at Re No. = 1.153 
with 1 Hz Imposed Unsteadiness 
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Figure 94. 
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RMS Velocity Profile at Re No. = 1826 
with 1 Hr Imposed Unsteadiness 
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RMS Velocity Profile at Re Wo. = 2026 
with 1 Hz Imposed Unsteadiness 
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RMS Velocity Profile at Re No. = 2086 
with 1 Hx Imposed Unsteadiness 
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Figure 142. 
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Figure 143. 
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Figure 144. 
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Figure 145. 



178 



PHASE ANGLE (doa 



RUNP8BR 93089.2106 



o 

CD 



w 

\ 

£ 

0) 

ID 

m 



J 

Ld 

cn 

r 

a 



w 

\ 

\ 

£ 

cn 

cn 



hi 

> 



e» 

> 

d 



0 


t- 

t 

i 

u~ 

L 


in n- 


i 


0J CD 


» 

J- 


in ru 


h 


0J 


c 

» 


0 


• # 


r 

r* 


■ 


4 - 


# > 


r- 


u 


j- 

u 


• a: 


a 




f~ 


f- 


tn 


i-* 

l 


w 


i 


\ 


i 

f- 


N £ 




I CD 


f* 






• 


r 


• CD 


r 


0 


f 


LI • 


i- 

i — 


a: _J 


r 


Li. Id 


> 


r 


• 


♦- 


<J Y 


r 


cn _i 




O ~J 


*“ 


CQ 


r' 

1 


146. 


~p-r t r 



ttjtt rTjTTrrp rrrrir n [ rr r rj ^ 



1 

1m- 

-~toj 

-jf*3 

-1 

Igd 

-103 

1™ 



c\j 

— LD 
-4pJ 

j 

1 

I 

-I 

ItD 

-*(M 



1 

J§ 



J. 

-j 



"~CD 

-*C3 



_i™ 

J 



CD 

cn 



ujaiu.Li.txU.iiJ.) 1 j iu 1 u.iUa.Lat-liijaJ iU-Lljaixifo 

3I % 21 • 80* MT 0 

TJA -3AU/lbH. n SWd 



/■N 



179 



PHASE ANGLE (deg 



Peak to Peak Magnitude of Phase Average 
RMS Velocity as a Function of y/d 




• 7 - 

O 



CO 

d 



r- 

r~ ; 



*u 

O 



L~! 

O 



K« 






Figure 147. 
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Figure 148. Re - 1136 
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Figure 149. Re = 1153 
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Figure 150. Re = 1278 
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Figure 151. Re - 1328 
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Figure 152 . Re = 1395 
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Figure 153. Re - 1395 
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Figure 154. Re = 1428 
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Figure 155. Re = 1428 
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Figure 156. Re = 1586 
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Figure 157. Re = 1586 
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Figure 158. Time Sequence at Re = 1428 
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